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Abstract 

Quantum quench dynamics is considered in a one dimensional unitary matrix model with 
a single trace potential. This model is integrable and has been studied in the context of non- 
critical string theory. We find dynamical phase transitions, and study the role of the quantum 
critical point. In course of the time evolutions, we find evidence of selective equilibration for a 
certain class of observables. The equilibrium is governed by the Generalized Gibbs Ensemble 
(GGE) and differs from the standard Gibbs ensemble. We compute the production of entropy 
which is 0{N) for large N matrices. An important feature of the equilibration is the appearance 
of an energy cascade, reminiscent of the Richardson cascade in turbulence, where we find fiow of 
energy from initial long wavelength modes to progressively shorter wavelength excitations. We 
discuss possible implication of the equilibration and of GGE in string theories and higher spin 
theories. In another related study, we compute time evolutions in a double trace unitary matrix 
model, which arises as an effective theory of D2 branes in IIA string theory. We find similar 
equilibrations and dynamical transitions in this matrix model. The dynamical transitions are 
related to Gregory-Lafiamme transitions in string theory and are potentially connected with 
the issue of appearance of naked singularities. 
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1 Introduction and Summary 

There has been a lot of progress recently in understanding dynamical processes in statistical 
mechanics models as well as in large N gauge theories. An important development in statistical 
mechanics has been an extensive study of equilibration in integrable systems pQ. In these studies 
it is found that equilibration happens only for a certain class of observables (whose characteri- 
zation remains an open problem), and when it does, it is not described by the standard Gibbs 
ensemble but rather by a generalized Gibbs ensemble (GGE) which has 'chemical potentials' for 
each conserved charge. There are a plethora of examples in string theory in which integrable 
field theories appear, including many which have gravity duals (see below). A question natu- 
rally appears whether we can understand equilibration in these models via GGE, and if so, how 
to interpret these in gravityll] One of the purposes of this paper is to initiate a study in this 
direction by studying Matrix Quantum Mechanics (MQM) models, which appear in a variety 
of contexts in string theory, as discussed below. 

Besides the above issue, MQM offers simple toy models to study dynamical phase transi- 
tions. Compared to some of the previously studied statistical mechanics models, MQM has an 
advantage that it has a simple large N limit {N is the rank of the matrix) which is described 
by a semiclassical Fermi liquid. The dynamics of the eigenvalue density of the matrix (whose 
gaps characterize the thermodynamic phases) is inferred in this limit by the classical motion 
of these droplets. In this paper, we study dynamics of a quantum quench, both for finite N 
(numerically) and in the semiclassical large N limit. 

The MQM (matrix quantum mechanics) we will study is the following % 



Z = I DU{t) e*^^, S = Idt 



lTT{m\')-v{u) 



(1) 



^In AdS/CFT and similar studies, thermal equilibrium has often been described in terms of black holes, 
which normally possess only a few parameters. A GGE, with an infinite number of chemical potentials, presents 
a novelty vis-a-vis this standard paradigm (see Section 2] for a detailed discussion). 

^Wc will restrict ourselves to excitations in the singlet sector of the model. This can be precisely achieved [5] 
by coupling to a gauge field Aq: dtU — >■ DtU, Dt = dt — iAo or by restricting to low energy dynamics [3] which 
cuts off the non-singlets. For our purposes, we may consider (fT9|) and (p4|) as the definition of the single-trace 
and double-trace models respectively. 
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We will focus on two cases: 

(i) the single-trace model, given by V{U) = - (TtU + Trf/^) , (2) 

and (ii) the double-trace model, given by V{U) = — — {TiU){TrU^). (3) 

Here U{t) is a unitary matrix of rank A^. The model (j2]) was first introduced in [2] in the 
context of 2+1 dimensional gauge theory. We will also discuss a related hermitian matrix 
model (l24l) which shares similar universal features with ([2]). We will, discuss in Section HI string 
theory/gravity duals of these models, in particular the non-critical 2D string theory and the 
near-horizon D2 brane geometrj^. 

We list below a few highlights of our result^: 

(i) We find evidence of equilibration [f| in the time-evolution of certain observables {e.g 
moments of the eigenvalue distribution) in both the MQM examples above. In case of the single- 
trace models, which is integrable, we find excellent agreement of the equilibrium configuration 
with that predicted by GGE (see Fig. [1]). We also give simple examples of observables which do 
not equilibrate. The observables that equilibrate are found to be the natural ones in the string 
theory/gravity duals, thus pointing to a possible characterization of them using gauge/gravity 
duality. 

(ii) We are able to give an interpretation of the equilibration in a special case where an 
initial long wavelength perturbation dissipates into progressively shorter and shorter wavelength 
perturbations, as in the Richardson- Kolmogor off cascade in turbulence. (See Fig. [3]). 

(iii) A surprising conclusion of our semiclassical analysis of the quantum quench is that 
dynamical transitions from the gapless phase to a gapped phase is possible, but not the other 
way around. In terms of the semiclassical Fermi liquid description (see Section 12.4.11) the 
interpretation of this is that Fermi liquid droplets cannot split. In the gravity context, such as 
in a Gregory-Laflamme transition between a black string and a black hole, these observations 
have a potential relation to a change of topology of the horizon. See Section 1431 for more details. 

The rest of the paper is organized as follows. In Section [2] we discuss the integrable single 
trace model and its dynamics, including equilibration and GGE, entropy production, energy 
cascades and dynamical phase transitions. In Section [3] we discuss the double trace model, 
including its relation to gauge theory, dynamical phase transitions and equilibration/oscillations 
for various initial conditions. In Section |4] we discuss implication of our results for the string 
theory /gravity duals of these models and possible generalization to a larger class of models. In 



^For other, recent, discussions of quantum quench in AdS/CFT, see [4] and references therein. 
^ The movies for our time evolutions are available on 



http : //www2 . yukawa . kyoto-u . ac . jp/-mtakeshi/MQM/ index . html 

^For some other works on thermalization in matrix models, see Refs [SI [6l [7], which consider BFSS and 
pp-wavc matrix models. 
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Appendix |X] we discuss the fermion picture of the matrix models, and the semiclassical Fermi 
hquid description at large A^. In Appendix [B] we discuss the notion of the generalized Gibbs 
ensemble in integrable models and its application to the single trace matrix quantum mechanics. 
In Appendix O we discuss the formalism of quantum quench dynamics in these models, and 
present explicit formulae for computing various relevant time-dependent quantities. In Appendix 
[D]we discuss a common limiting case of both the single-trace and double-trace models in which 
the fermions are free and are without a potential; the density profile has a nontrivial nonlinear 
dynamics and shows, in this case, a power law decay to equilibrium. In Appendix [E] we present 
some details of our numerical analysis used to compute time evolution. 



2 Single trace model 

The single trace model ([2]), as explained in Appendix |X] and in the references cited there, can be 
equivalently formulated in terms of noninteracting fermions moving in one spatial dimension 
in an external potential V{6) = a cos 6*. The hamiltonian, in the second quantized form, is (see 

(USD) 



2A^2 



H = jde ^\d,t) 

pio, t) = ^i^\e, t)^{e, t), jde pie, t) = i. (4) 

The fermion density p{e,t) is the same as the eigenvalue density of U{t). Depending on the 
parameter a, there are two different static phases (ground states) in the system: (i) gapless for 
a < ac = n^/GA, and (ii) gapped for a > ac, where the 'gap' refers to that in the equilibrium 



value of p{e) (see Fig. H] 



or various shapes of p{e)). The quantum critical transition at Oc is a 



third-order transition |2] ^ (see below ( 12^ 1 

In what follows, we will analyze various types of dynamics of the above system: 

(i) We study time evolution of some arbitrarily chosen initial state (corresponding to a suitable 
density profile) under the hamiltonian (jl]). 

(ii) Quantum Quench Dynamics (QQD) (see Section O) for details): we start, in the far past, 
with the ground state of the Hamiltonian (jl]) with an initial parameter a = af, we then change 
a suddenly to aj at t = 0, and study the t > evolution of this initial state and various 
expectation values under the hamiltonian (jlj) with parameter a = af. We can regard this as a 
special case of (i) where the initial state is provided by the ground state of the system in the 
far past. 



^This is similar to the Gross- Witten-Wadia transition [5J H] of zero-dimensional single- trace unitary matrix 
model. 
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2.1 Integrability and equilibration 

Since the hamiltonian (jlj) describes free fermions, it is clearly integrable. Indeed we can find the 
conserved quantities as follows (see Appendix [B] for details). Let us expand the fermion field as 

oo 

^(^) = X] Cm^miQ)-, hip^{0) = em^rn, (5) 
m=l 

where ipm is the m-th eigenfunction of the single particle hamiltonian h ffTSjl with the eigenvalue 
em and Cm is the corresponding annihilation operator. The fermion occupation number operators 

Nm = ci,Cyn, (m = l,---), (6) 

are clearly all conserved, and all independent in the — oo limit. An equivalent basis of 
conserved charges is the set Ip defined in ( l36l) in Appendix [Bl For finite A^, the set {/i, ...,In} 
provide an independent set of conserved charges. 

A first guess about a free system such as this would be that no observable can possibly 
equilibrate since there presumably cannot be any dissipation. Indeed, such a guess would seem 
to have a lot of merit. Firstly, there are infinite number of charges, such as the N^, which, given 
any initial state however complicated, clearly stay constant as the system evolves. Secondly, 
it is easy to construct a whole series of observables which do not remain constant but keeps 
oscillating for ever; as an example consider the operator 

Om,n{^) = cl^Cnexp[-i{em - e„)t//i], m^n. (7) 

Clearly, given any state the expectation value of the above operator evolves as (using 

the Heisenberg picture) 

(6^,„(t)) = (*(0)|6^,„,(t)|*(0)) = K»,nexpH(e„ - e^)t/h], Vm,n = (^(0)|4c„|*(0)). (8) 

Thus, the expectation value keeps oscillating ad infintum (unless Vm,n = in the given initial 
state, in which case the value remains zero. It is obvious how to generalize such examples, e.g. 
by combining many creation-annihilation pairs. 



Behaviour of p{9, t): In the light of the above discussion, we find surprises when we consider 
the time evolution of the density p{0,t) or its various 'moments', defined by 

{Pn{t)) = j de{p{e,t)) cosine). (9) 

In (l56l) a finite- A^ expression for {p{6,t)) is presented for a quantum quench (defined briefly 
below (jl]) and discussed in detail in Section [C]). Using this expression, it is easy to compute the 
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time-dependence of various moments pn{t)- In Fig. [T] (a) we have plotted pi{t) for = 120. 
In the same figure we have also plotted the — > oo limit of pi{t), computed using (|9]) and the 
expression fl32|) . the latter being derived using the large semiclassical limit of the fermion 
theory (j4j). From these plots, the evidence of equilibration is quite clear; while for finite A^, 
there appears to be relaxation followed by a recurrence (with the recurrence time increasing 
with A^ as shown in Fig. [1] (b)), at large A^, pi{t) appears to relax to a certain value. 




20 40 60 80 100 120 140 o 100 200 300 



(a) pi(t) and GGE vs. Gibbs ensemble (b) Poincare recurrence 

Figure 1: The time evolution of the first moment \{pi{t))\ (defined in ^) is shown. We take = O.Suc 
initiahy and change it to a/ = 1.2ac instantaneously at t = 0. (a) The (red) dashed line shows the result for 
N = 120, which appears to show Poincare cycles. The Poincare cycles seem to grow linearly with TV as shown in 
(b). The (blue) solid line in (a) is for iV = oo (computed using the droplet formalism, see Section 1X1 and IE)) : it 
shows equilibration and absence of Poincare cycles. The mean value of \{pi{t))\ for finite TV and its asymptotic 
value for N ~ oo are seen to have excellent agreement with each other and, in turn, with |(pi)gge| as computed 
in GGE, (see Section [B] for details of the computation), whereas these values all differ significantly from the 
ensemble average in the standard Gibbs ensemble. We have verified that the N dependence of the (yOi)GGE is 
small. 

Understanding of the equilibration in terms of GGE: Since in view of the discussion 
around Eqn. ([H]), the above equilibration is surprising, it is useful to look for similar phenomena 
in other examples. Fortunately, in recent years several examples of selective equilibration have 
been found in integrable systems (see [1] for a review), in which a certain class of observables 
has been shown to equilibrate in such integrable systems (which include the hard core bosonic 
lattice and the transverse field Ising model), and the equilibrium configuration is characterized 
not in terms of the standard thermal (Gibbs) ensemble, but by a Generalized Gibbs Ensemble 
(GGE) which keeps track of the infinite number of conserved quantities by means of an infinite 
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number of chemical potentials! In our present model, the GGE is defined by the density matrix 

0GGE = -7^— exp ( - V /i^iV„ 

where /i^ is the chemical potential corresponding to the conserved fermion number ^ 
(see sections [B] and [C.ll for details). Using this density matrix, we are able to compute the 
postulated equilibrium value of the density (p(6'))gge = Tr (fPcGEPl^*)) and hence, through ([9]), 
that of the various moments (pn)GGE- In Fig. [T] we have shown the value of (pi)gge; we find that 
the large t asymptotic value of {pi{t)) computed at large fits very well with the GGE- value. 
For comparison, we have also displayed in the same figure the value of pi computed according 
to the Gibbs ensemble, which differs significantly from both the asymptotic value of {pi{t)) as 
well as (pi)gge- We have verified these statements for higher moments P2, Ps etc. as well. We 
find, therefore, that the moments p„(t) belong to the class of observables in our system which 
and the equilibrium configuration can be quantitatively understood in terms a 
GGE. It is also interesting to note that there are many observables such as ([7]) which never 
equilibrate. We emphasize that the density variable p{6, t) which does equilibrate and satisfies 
the GGE hypothesis (144|) . has a direct interpretation in terms of the string theory dual (at least 
in case of the hermitian model f l2^ ). 

2.2 Entropy 

The phenomenon of equilibration is intimately tied to the idea of coarse-graining and entropy 
production. To understand this, let us rewrite the GGE hypothesis (Hil) as follows: 

Tr(ftrue 0{t)) ^ Tr(f?GGE 0(t)), as t ^ 00, (10) 

where ^tmc = |^(0))(^&(0)| is the density matrix corresponding to the initial state at t = 0. 
Note that the von Neumann entropy S'truc = — Tr Qtmc In ftruc vanishes since Qtme corresponds 
to a pure state. On the other hand, the von Neumann entropy ^gge = — Tr f?GGE In f?GGE, 
computed in ( l43|) . is non-zero. Thus, clearly, ( fTOl) cannot be interpreted to mean Qtrucit) = 
e~*^*f?truoe*'^* *— >° Qgge ! The correct interpretation of ( !T0|) is that the class of operators 
0{t) are suitably coarse-grained to smooth out the difference between Qtme and ^?gge after a 
sufficiently long time. From the discussion in the previous subsection, we see that the density 
operator p{9, t) is one such operator. 



^The infinite number of chemical potentials can sometimes be thought of as a separate 'temperature' for 
every mode, see, e.g. [10]. 

*The present example provides evidence for relaxation for a quantum quench dynamics. However, we will 
find later the same phenomenon with other initial conditions; in particular see Section [D] where we provide an 
analytic proof of such relaxation in a special case. 
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Amount of entropy production: Using fl43l) and (135]) we can compute tlie entropy S'gge 
cliaracterizing tlie GGE. In Fig. [21 left panel, we show the result of this computation as a 
function of A^. From the plot we can see that S'gge oc This could be roughly expected from 
the fact that there are non-interacting particles in the system. In Section H] we discuss this 
0{N) entropy in terms of various string theory duals. 



0.30 Entropy /A^ 
0.25 




(b) tti > ttc 



Fi gure 2: Sqq^/N vs ay. These results indicate that S'gge scales with N. We can also see that they are 
proportional to \ai — af\. 



2.3 Energy Cascade: an interpretation of the dissipation in the in- 
tegrable model 

To get more insight into the phenomenon of equilibration, let us consider the special case of 
V = Q (i.e. a = 0) in (j4]). In stead of a quantum quench type initial configuration, let us in this 
case start with an initial state |\l/o) which corresponds to a sinusoidal semiclassical configuration 



see 



02B|) for notation^ 



p(^,t = 0) = — (l + 26cos(n^)), V{d,t = Q) = 0. (11) 
27r 

The initial condition implies that only p„ is excited at t = 0. By using the technique detailed 
in Appendix ini we can compute the time evolution of the system (see Figure [T3|) . and compute 
expressions for the moments Pn{t) analytically. For instance, the late time behaviour of the first 
two moments (for the choice n = 1 in ( ITTl) ) is given by 

PiW ^ T^\/4:Cos ^ ) cos ( ) , (t^oo), (12) 



^We have also verified that in units of (|4]), the conserved energy of the system is i? = 0{N) and the 
temperature is T = 0(1). 

"'^'^This kind of initial condition is chosen here only for simplicity; we expect the same qualitative behaviour 
for the quantum quench initial conditions as well. 
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and 




Figure 3: Time evolution of moments of the eigenvalue density for the free model (sec Appendix ID] for details). 
Note the time of appearance of the first peaks. Lower frequency modes get excited and die out first. Higher 
frequency modes get excited and die out later. There is a flow of energy from longer wavelength modes to shorter 
wavelength modes, as in the Richardson-Kolmogoroff cascade in turbulent fluids. 



It would appear surprising at first sight that ALL the moments die away!, leaving us ap- 
parently with the question 'where does the energy go?' Of course, since we have a conservative 
model, even integrable to boot, energy cannot really 'go' anywhere. It turns out, just as studies 
in turbulence and various other areas have taught us, that there is a flow of energy from one set 
of 'modes' to another set of 'modes'. To be precise, as the longer wavelength modes die, shorter 
wavelength modes start getting excited, and as they too die away, even shorter wavelengths 
get excited, etc. This phenomenon is similar to the Richardson-Kolmogoroff type cascade in 
turbulence, and provides a mechanism of dissipation in our integrable system. 

2.4 Dynamical phase transitions 

We mentioned above that if we tune the parameter a in (jlj), say from zero upwards, we encounter 
a third order quantum phase transition at a = = 7r^/64 from a gapless phase to a gapped 

^^The power-law relaxation may be a special feature of the simple a = model. In the more general case, 
such as in Fig. [1] the relaxation appears to have a characteristic time scale. 
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phase. What happens if we make this change of a dynamically, say from < to > ac? If 
we prepared the system in the ground state initially (corresponding to a = Oj), does it eventually 
relax to the final ground state (corresponding to a = a/)? Do we observe dynamical appearance 
of a gap at a certain time? What role does the critical point Oc play, if any? What happens if 
we make the change from at > ac to aj < acl 

The general setup for addressing these questions is to consider a dynamical variation with a 
tunable 'ramp speed', e.g. in the Kibble-Zurek transition (see, e.g., the review [1]). However, in 
the present work, we will consider two special cases, namely (a) the quantum quench dynamicJ^ 
discussed in the previous section, in which the change from ai to a f happens instantaneously, 
and (b) the adiabatic case, in which the change happens at a rate slower than the scale set by 
energy level spacings. We will first address the question of appearance/disappearance of a gap 
in the final state in case we cross the critical value ac from an initially gapless/gapped phase. 

2.4.1 One-way phase transitions 

Fig. m shows, for specific values of ai,af, the large t asymptotic density p{0,t) for quantum 
quench dynamics (QQD), calculated from fl56|) . as well as the same quantity for adiabatic 
dynamics (where the time-dependent wavefunction is approximated by the instantaneous ground 
state), as calculated in ( l6T]) . We find that there is a rather remarkable 'time-asymmetry' in the 
QQD, viz. that while a gapped — )■ gapless transition is possible, a gapless — )■ gapped transition is 
impossible. In other words, in the context of QQD, while a gap can close dynamically, it cannot 
open dynamically!^ 

In Fig.[5]we present more details about the gap at varying values of the a-parameter. As seen 
in Fig. m for a > ac the gap in the ground state density p{9) opens around 9 = (because the 
potential a cos 6* draws the fermions from 6 = towards ^ = vr, as explained below fl23|) ). Thus, 
the value of the density at ^ = can be regarded as an order parameter which is non-zero in 
the gapless phase and zero in the gapped phasJ^. In Fig. |5l we present this order parameter for 
various values of aj for QQD after relaxation and compare it with the ground state. The one- 
way nature of transition is quite clear in this plot as well. Note also the progressive departure 
from adiabaticity as a/ is varied, starting from a,, across ac to the other phase (this is more 
prominent in the left panel where we study the appearance of a gap). We will have more to say 



^'^ In this section, we use GGE to evaluate the asymptotic values of the observablcs in the QQD, although, at 
finite N, the actual values keep oscillating slightly around the GGE predictions as shown in Fig.[T](b). 

^■^Of course, this is not a real time-asymmetry in the sense that, for < Oc < o/ if we start from the rather 
complicated, gapless, final pure state, and reverse the direction of time, we will indeed get back the initial gapped 
pure state. However, in the context of QQD, we always start from a ground state, in either direction, leading to 
the apparent time-asymmetry. 

^^These statements refer to the large N limit. For finite N, or in a double scaled limit, a non-zero fermion 
density exists in the 'gapped region'. 
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-3-2-1 1 2 3 0-3-2-1 1 2 3 



(a) tti < ac < af (b) ai > > af 

Figure 4: A plot of density {p{0)) vs 9 at N = 120. In the left panel, a/ > Oc > a^; the green curve shows 
the initial gapless density as a function of 9, the blue curve shows the gapped density corresponding to the final 
ground state for a = af, and the red curve shows the large t asymptotic value of p{9,t) in quantum quench 
dynamics. In the right panel, Uf < Qc < a,i, the initial density is gapped, whereas both the final configuration in 
the actual dynamics as well as the final ground state show a gapless phase, showing the presence of a dynamical 
phase transition in this case. 



shortly about the departure from adiabaticity. 



P{0 = 0) p{9 = 0) 




(a) tti < ttc (b) tti > ttc 

Figure 5: A plot of p{9 = 0) for various values of a/ at = 120. In the left panel, we begin in the gapless phase 
with aj = ai < Oc, and consider a continuous sequence of QQD evolutions in which a/ is progressively increased 
across criticality to a/ > Oc : the 'adiabatic' curve, calculated according to (j6ip . shows transition to a gapped 
phase (with p{9 = 0) ^ 0), whereas the actual dynamical evolution under quantum quench dynamics (QQD) 
shows that no gap appears dynamically {p(9 = 0) remains non-zero). In the right panel, both the adiabatic and 
QQD curves show disappearance of the initial gap, confirming the presence of a dynamical gapped — > gapless 
transition. Note that the gap does not disappear at Oc even in the ground state (adiabatic case). This is due to 
a finite N effect. 
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2.4.2 A topological interpretation of the one-way transition from Fermi liquid 
drops 

The apparent irreversibility in dynamical phase transitions, found above, turns out to have a 
beautiful semiclassical interpretation from the semiclassical Fermi liquid picture. For simplicity, 
let us consider the large limit of the single-trace hermitian model (124|) . which corresponds to 
replacing the potential a cos 6* in (jlj) by —a(P'j1 (defined together with a hard cut-off Q = ±6*^)- 
Since in these models h = 1/N, the large N limit, given by fl25]) . is given by the semiclassical 
picture of a Fermi liquid in phase space, which is described by a phase space density u{9,p,t) 
which is either or 1, and the filled regions u = 1 are described as 'droplets'. There is a simple 
description of quantum quench dynamics in this picture, as already alluded to in Section 12.11 
and described in detail in Section \M 

As we describe in Fig. [6l the time evolution from the gapless phase towards the gapped 
phase involves starting with a single connected Fermi droplet in the Fermi sea. The evolution of 
the droplet follows the equation f l30p . It is not difficult to see that the hyperbolic motion of the 
fermions (in which fermions belonging to different hyperbolas move at different speeds) leads to 
squeezing and stretching of the droplet, although it can never split. Since p{6, t) = J ^u{6,p, t), 
which can be interpreted as a kind of projection on to the 6'-axis, the impossibility of droplet 
splitting in phase space translates to the impossibility of a gap-opening transition in ^-space. In 
other words, a continuous u{6,p,t) means that there will always be some fermions in the phase 
space along the p-axis, hence p{6 = 0) will remain non-zero. We explicitly see this in part 1 of 
the figure. On the other hand, when we consider time evolution from a gapped phase towards a 
gapless phase, even though the initially disconnected Fermi sea remains disconnected in phase 
space (again working with the same evolution equation fl30l) ). we find that the projection on the 
6'-axis can become gapless, thus ensuring the presence of a gapped — )■ gapless transition. 

Relation to Gregory-Laflamme transition: Note that the above argument is robust and 
ensures that the asymmetric nature of dynamical transitions will not depend on the details 
of the potential V{U) in ([l])(for example, it persists for the double trace potential, as we will 
see in Section [3l). Now the choice of V{U) dictates the rate at which the phase space fiuid is 
squeezed due to the dynamics of individual particles; what happens when V{U) is such that this 
squeezing effect is very strong? In this case, the neck will be stretched thin and its thickness will 
shortly become 0{1/N) = 0{h). At this stage, the semiclassical droplet description will clearly 
be violated, and we will need to use the full quantum mechanical description. As discussed in 
Section 14.51 a gap-opening transition in the double trace matrix model is related to a Gregory- 
Laffamme transition in the dual string theory. In a GL transition, classical general relativity is 
not valid when the size of the 'neck' (of a black string or of a solitonic string) becomes Planck size. 
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In case of the black string — )■ black hole transition, the breaking of the neck involves appearance 
of a naked singularity. In view of the discussion above, this singularity would be related to the 
size 0{1/N) neck region in the matrix model. Thus we expect that the understanding of matrix 
models would be important to understand the issue of the Gregory- Laflamme transition. We 
will discuss this issue further in Section 14.51 




(I) ai < ttc, af > ac. 




(II) tti > ac, ttf < ttc- 

Figure 6: Quantum quench dynamics of the filled Fermi sea in the hermitian model ([M)) . The dashed lines are 
p = ±y/af9, which are the separatrices for the single particle motion. In (I), the left panel shows the shape of 
the initial filled Fermi sea, corresponding to the ground state for a = ai < ac- This is a gapless phase, hence the 
Fermi level is above the top and the droplet is a single connected one. The right panel shows an intermediate 
time-evolved configuration. After the quench, the modified hamiltonian stretches the droplet, but it cannot split 
it, hence the eigenvalue density, which is a certain projection on the horizontal axis (see (|26p ). also cannot split. 
In (II), left panel, the initial a ~ ai > ac is in the gapped phase. The Fermi level is below the top, hence the 
left and right 'wells' are both filled up to the Fermi level, while remaining disconnected. Hence the filled Fermi 
sea has two disconnected segments (therefore, the eigenvalue distribution is gapped). The right panel shows the 
time-evolved configuration (according to a = ay < Uc), which is such that although the droplet still consists 
of two disconnected parts in phase space, the projection on the horizontal axis is gapless, hence the eigenvalue 
density is gapless. 

2.4.3 More on departure from adiabaticity 

In this subsection, we come back to the issue of the role played by the quantum critical point in 
quantum quench dynamics (QQD). This has been discussed in several contexts: e.g. (a) how the 
asymptotic density of excitations scales as a function of the parametric distance from criticality 
when the initial hamiltonian is at criticality (see [H]), or (b) how the asymptotic QQD value 
of some observable or order parameter differs from the adiabatic value, as a function of (in 
our notation) af — (for given fixed Oj). In the above subsections, we found many examples 



13 



of differences between QQD and adiabatic evolution. We explore this aspect more in this 
subsection by studying the behaviour of pi. We find that as aj is taken further away from Oj 
towards the other side of the quantum critical point, the departure from adiabaticity becomes 
appreciable from near aj ~ Oc- We have found similar behaviour for the higher moments P2,P3 
etc. as well. This phenomenon of departure from adiabaticity around the critical point is similar 
to that reported in [12] in the context of QQD in a Hubbard model. Indeed our Fig. [7] is quite 
similar to Fig. 3 of [12],^ This points to some possible universal role played by the quantum 
critical point even in dynamical phenomena. 
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(a) flj < ttc 



(b) tti > ttc 



Figure 7: A plot of (the asymptotic value of) pi as a function of a/. The setup is the same as that of Fig. [5] In 
the left panel, we start with the gapless phase aj ~ ai ~ 0.8ac, and consider QQD experiments for progressively 
higher values of a/ till we go across af = Qc towards a/ = 1.2ac- The blue curve shows the adiabatic value of pi 
(which corresponds to the ground state of the a/-hamiltonian), while the red curve shows the large t asymptotic 
value of pi{t), which agrees with the GGE result. The two curves start showing significant departure from 
around af = Uc- In the right panel, we start with the gapped phase with a/ = a.i = \.2ac and end up with 
af = O.Sflc; the colour coding is similar. Once again we find that departure from adiabaticity starts from around 
af = flc. This figure is similar to, e.g.. Fig. 3 of |12] . where a similar pattern of departure from adiabaticity is 
reported for quantum quench in a Hubbard model. 



3 Double trace model 

Like the single trace model, the double trace model ([3]), can also be described (see (134|) ) in terms 
of fermions moving in one spatial dimension. The hamiltonian is given by 



H = Jde l^deij\e)de^{e,t) + Hi^u 



H.nt = j;^ I dtdd d9'[i;\9, t)^ie, t) cos{e - e')^\9\ t)^{e\ t)]. (14) 



^^The departure from the adiabatic curve does not exactly coincide with criticality in [12] or in our work here, 
an effect which is more pronounced for higher moments. 
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The difference from the single trace case is that the fermions are now interacting with a mutual 
^ cos(6' — 6') potential. Since this is an interacting system, it is not obviously integrable; we will 
not explore the question of integrability in this work. 

Before going to explore dynamics in this system, we will briefly review its connection to to 
2D adjoint scalar QCD and string theory. 



The 2D adjoint scalar QCD: Consider the following two dimensional gauge theory on a 



S 



dt I dx Tr 




1=1 i,j J 



(15) 



Here L is the period of the S^. This theory is related to a dimensional reduction of D + 



2 dimensional pure Yang-Mills theory (w 
compactified on a Scherk-Schwarz circle^^^ 



rich, for D = 8, can be related to D2 branes 
The phase structure of this model is as shown in 
Fig. IHl At sufficiently large /3, (at least above the critical line CBDO), the adjoint scalars would 
have a mass gap A and are integrated out through the 1/D expansion [13] (see also [531 E])- 
The theory is then described by the following effective action 



S/DN^ 



dt 



71=1 



Ki(nAL) 



n 



— Trf/" 
N 



(16) 



where U{t) = P exp[i f Ai{t)dx], and Ki is the modified Bessel function of the second kind. 
In the following subsections, we would mainly consider dynamical phase transitions near the 
critical line AB in Fig.|8l In this case, we can ignore the higher modes {n > 1) in the infinite sum 
in the potential term, and obtain the double-trace model ([3]). Here ^ in ([3]) or ( iMj) is a function 
of the parameters in the QCD, and importantly it monotonically decreases as L increases [13] . 



Phases: The parameter ^ determines the phase structure of the theory. For high enough ^ 
(^ > ^2), the attractive potential Smt dominates and the system is gapped. For low enough ^ 
(^ < ^2), the kinetic term dominates and the system is gapless. See Fig. [8] (b) for more details. 

3.1 Dynamical evolution 

In case of the double trace model, several solutions have been derived for a given C, (> ^0) 
as shown in Fig. [8] (b) but the ground state is always unique. Hence we naively expect that 

^^If we derive the model (|15p from the higher dimensional gauge theory or D branes, additional terms are 
induced through quantum loop effects. For this reason, the coupling at the two dimensional gauge field g and 
the commutator interaction g' have been distinguished in the model (|T5|) . The most relevant term in these loop 
corrections in weak coupling is the mass term for the adjoint scalar. We omit this mass in this paper as argued 
in [ig. 
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TtU yi 




^ Phase II (gapped) 

^ Phase IV (non-uniform) 
Tr(7 ^ 




Phase I (uniform) 
TrC/ = 

Phase III I 
TrC/ 7^ 



(a) Phase structure of 2D adjoint scalar QCD 



(b) A plot of free energy around line AB of 
diagram (a). 



Figure 8: (a) Phases of the 2D QCD model (fT5|) are shown for various vahies of the inverse temperature /3 
and spatial size L. The phases are characterised by the expectation values of the temporal Polyakov loop V and 
the spatial Wilson loop U . We are interested in the limit of zero temperature (non-compact time direction), 
in which the Polyakov loop 1/ = 0; the only phase boundary of interest here is the line AB^ across which the 
spatial Wilson line U changes. The theory in this region is described by ([3]) or (fH)) . (b) The energy H ^ E oi 
(fT4|) . computed using large TV Fermi liquid picture, is plotted for various values of ^ (see [13] for details). To 
connect with diagram (a), note that ^ is a monotonically decreasing function of L. Solid lines represent the 
stable phase (the solution with the lowest value of E for a given ^) as well as metastable phases, while the dotted 
lines represent unstable phases. At ^2 = 0.237 there is a Ist-order transition from a stable uniform distribution 
to a stable gapped one. Below ^0 = 0.227 the metastable gapped phase does not exist, similarly above ^3 = 1/4 
the metastable uniform phase does not exist. At ^1 = 0.231 there is a Gross- Witten-Wadia (GWW) transition 
among two unstable phases. 



the unstable or meta-stable solutions may evolve toward the ground state if we add sufficient 
perturbations to these solutions. We will explore the question of possible dynamical phase 
transition as well as the question of relaxation in this process. In this section, we use the 
semi classical approximation. (See the details in Appendix |El) 

3.1.1 Dynamical evolution from the uniform gapless phase 

We will consider the time evolution of the uniform phase. In terms of the droplet picture 
(Section El), the uniform phase is given by 

p(^) = l/(27r), V = 0. (17) 

This configuration is always a solution of the model (fT4|) . since the force term vanishes for 
particles located uniformly on the circle. In order to have non-trivial dynamics, we add a small 
perturbation to this configuration. 

Before we discuss the actual results of the evolution below, let us pause and ask what to 
expect. If ^ in f|T^ was small enough, say well below ^2 = 0.237, for such values of C, the uniform 
phase is the ground state. Thus, we expect that for small enough amplitude of the perturbation, 
the system will revert back to the uniform phase. Indeed for ^ = 0, discussed in Section [2^31 we 
find this to be explicitly true. 
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On the other hand, we are interested here in values of ^ above .^2 where the ground state is 
the gapped one. Thus, we will expect the perturbation to be unstable towards formation of a 
gapped phase. 



Order parameter : Since all moments ([9]) of the density vanish in the uniform phase, any of these, 
say pi, can be considered as an order parameter, which would change from pi = — pi 7^ 
if the distribution changes from uniform to a non-uniform solution. Note that we consider Z2 
symmetric {6 ^ —6) perturbations only in this article. Then we can take pi = TiU/N by a 
gauge choice. 

We will now discuss the actual result of the time-evolution of a perturbation according to 
(1741) . In Fig. [9l the result is presented for ^ = 0.260, at which the uniform solution ( |T7|) is 
unstable. Then |pi| develops to a non-zero value (see the right panel of the figure), and we 
obtain a non-uniform distribution. However the distribution does not develop a gap, unlike the 
naive expectation above. Indeed, just as in Fig. El part (I), the gapless — gapped transition is 
once again not possible since the droplet in phase space does not split. In addition to all this, 
we also find evidence for equilibration in |pi|(see the right panel) which indicates relaxation of 
the phase space configuration to one with a thin stable neckF^ 



pi 




t = t = 200 pi(t) 

Figure 9: Time evolution of the uniform solution ([T7)) with a small perturbation at ^ = 0.260. The left panel 
shows the initial configuration. The configuration approaches a gapped phase, but develops a thin neck, like in 
Part (I) of Fig. inland the gap never develops. The right-most panels show the behaviour of the order parameter 
pi with time. It shows evidence of an equilibration. 



3.1.2 Time evolution from the gapped phase 

The discussion is analogous to the previous case. We consider a meta-stable gapped state at 
^ < ^2 (Phase III in Fig. |H1 (b)) and perturb the system. 

The result at ^ = 0.230 is shown in Fig. [TOl A dynamical gapped — )■ gapless transition takes 
place. Thus, like in the single trace case, we once again see a one-way dynamical phase transition, 

"'^'^The 'neck' avoids a gap-opening transition, and consequently avoids a dynamical singularity. It is important 
to explore this further, particularly in the context of the double scaled matrix model. For similar issues related 
to smoothening of singularities through 1/iV effects, please see [HI 115) . 
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gapped to gapless, but not gapless to gapped. This once again, has a clear interpretation in terms 
of phase space. 

Note that if the perturbation is not sufficient, the state reverts back to the original gapped 
configuration, since the configuration is meta-stable. 




Figure 10: The figure on the left shows a slightly perturbed meta-stable gapped distribution at t = 0. The 
value of ^ is 0.23. The figure in the middle shows a gapless distribution at t = 400. The figure on the right 
depiets pi{t) as it changes from 0.4 at i = to 0. 



3.1.3 Appearance of exotic states 

We have considered the time evolutions of unstable/meta-stable solutions by adding pertur- 
bations. However if we add sufficiently large perturbations, we observe appearance of exotic 
asymptotic states. See Fig. [HI e.g., where in the middle panel two blobs of fermions have 
appeared, which are moving towards the left and the right respectively. As time progresses, the 
two blobs keep executing a period motion and in coordinate space, this means that two density 
peaks approach, merge and separate out, and repeatedly pass through each other like colliding 
solitons (remember that in our model the fermions are moving in a circle). 



pi 




Figure 11: Time evolution of the unstable uniform solution ([TT]) with a large perturbation at ^ = 0.260. The 
partieles split into two well-defined blobs in phase space. The upper blob moves to the right while the lower blob 
moves to the left, and they keep moving past each other periodically. In coordinate space, the blobs overlap; 
this implies that two density peaks periodically pass through each other like colliding solitons. 



Another exotic time evolution can be seen if we consider the higher modes of the potential 
(fT6|) . We assume appropriate values of A and L such that ^ is large, consequently the gapped 
solution is stable and the uniform one ( IT71) is unstable. In this setup, the unstable uniform 
solution f lT7|) evolves as shown in Fig. [T2j In this case, the dynamical phase transition consists 
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of two steps. First the unstable uniform state decays to a 2-peak state (where the 'peaks' refer 
to those in the coordinate space density p{9)). This 2-peak state remains for some time, and 
finally decays to the one-peak state. We have also observed appearance of various multi-peak 
intermediate states if we change the initial perturbations. 

Indeed the existence of unstable multi-peak solutions for the small L region in Fig. [8] (a) 
and the possibihty of their appearance during the decays of the unstable uniform solution 
have been predicted in [16]. Our result confirms this conjecture. Interestingly similar time 
evolutions have been observed in the decay of black string in the Gregory-Lafiamme transition 
also [591 EI]- An unstable black string does not always decay directly to a single-black hole, 
which is thermodynamically most stable, but to multiple black holes in some cases. 



Figure 12: The appearance of the 2 peak state during the time evohition of the uniform configuration. We 
take ^ = 5.0 and AL = 1.0. The strong signal of p2 indicates the two peaks in the eigenvalue. First the uniform 
state decays to the 2-peak state around t = 1, and then it decays to one peak state around t = 15. (We do not 
plot the one peak solution here but it is similar to Fig. |9]). 

Note that the two-peak state that appears in Fig. [12] is qualitatively different from the 
breather-like asymptotic state in Fig. [11] The former state is approximately static (with two 
distinct, long-lived density peaks around 6 = ±7r/2) until it decays to a single-peak state, 
whereas the latter state is approximately stationary and the two peaks merge and separate 
out repeatedly in an oscillatory fashion. (See the movies for these evolutions from the link in 
footnoted) 

As we have seen in section [2] the system is integrable if the double trace interactions are 
turned off. One open question is whether this interaction breaks the integrability. The appear- 
ances of these exotic asymptotic/intermediate states may shed light on this issue. Especially 
the exotic asymptotic states in Fig. [TT] may imply an attractor structure in this system. We 
will come back to this issue in a future work. 
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4 Discussions 



The results of our paper have already been summarised in Section [T] and described in Sections 
[2] and [3l In this section, we will describe some related issues, speculations and outstanding 
questions. 

4.1 Duality to 2D string theory and W^o algebra 

The model defined by Eqs. (IMl) . for fixed a, is essentially the c = 1 matrix model. This 
model (in its double-scaled form: (a — ac)N = fi), has a string theory 'dual', namely the 
two-dimensional bosonic string theory, which is characterized by a massless scalar field (called 
'tachyon') on a fiat spacetime and linear dilaton background (see, e.g. the review of [H])- The 
tachyon field fiuctuation is given by a non-linear and non-local transformation of the matrix 
model scalar field fiuctuation (see, e.g., [191 1201 HI] for more details). The filled Fermi sea, 
p{6) = ^y^2{fi — V{9)), P = 0, is a solution of ([27]), and corresponds to the tachyon vacuum 
in string theory. In quantum quench dynamics, in the matrix model description, we start 
with an unstable density perturbation, which evolves in time under the restoring force towards 
the above-mentioned stable configuration. Since the motion is conservative (in fact, it has an 
infinite number of conserved quantities), the configuration cannot reach the stable configuration; 
in stead there is an oscillatory motion, in which lower frequency modes transfer energy to the 
higher frequency ones, leading to an energy cascade as discussed in Section 12. 3i^^ 

It would be interesting to map this oscillation to string theory tachyon modes by applying the 
transformation alluded to above. We expect that qualitatively it would be a similar phenomenon 
to the one found in the matrix model, with a finite energy extended configuration giving rise to 
proliferation of high frequency modes, asymptotically approximated by a general Gibbs ensemble 
(GGE) in spacetime. Note that in the 2D string theory too, there are an infinite number of 
conserved charges, which are the diagonal elements of an infinite dimensional algebra, viz. the 
Woo algebra (see [22]), thus making the GGE a natural choice of ensemble. It is an interesting 
question whether the GGE can be represented by a new field configuration such as a non-trivial 
solitonic configuration, or a black hole0 We discuss the black hole scenario in some detail in 
the subsection 14.31 below. 

In [23], a different, type 0, 2D string theory is proposed as dual to the (gauged) hermitian 
model (El]), where the singlet condition is ensured by a gauge covariant kinetic term Tr((9t + 
i[Ao,M])'^. In the context of the type OB theory, there are two massless scalars, the tachyon 

"'^^We show this in the unitary matrix model, but wc find similar behaviour also in the Hermitian matrix model 
which is directly identified as the c ~ 1 matrix model. 

^^It has been argued that inflaton oscillations around its vacuum configuration, which is analogous to the 
tachyon oscillations described above, can produce primordial black holes (see [23], for example). 
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and an axion. We expect many of the comments in the previous paragraph to go through in 
this case as welL 

4.2 Gravity duals of other integrable systems 

The phenomena of selective equihbration and emergence of GGE studied in this paper would 
be of interest in many other examples of integrable systems with gravity duals. We mention 
below a few of them. 

D1-D5 system: The D1-D5 system (reviewed in, e.g. [25]) is described, at a certain point 
P of its moduli space, by a free 2D CFT (with target space an orbifold of the type / S{N), 
where M = or K3, and = NiN^). The usual D1-D5 supergravity solution is located 
at a different point Q of the moduli space, which can be reached by introducing a marginal 
interaction A J O to the free CFT at P. There are certain nonrenormalization theorems which 
ensure the protection, under the A-deformation, of a number of quantities including black hole 
entropy and absorption cross- sect ion. In [26] a certain response function in the A = integrable 
theory is exactly computed and shown to agree with its gravity dual, namely the absorption 
cross-section of the two-charge black hole (at the A = Aq supergravity point). The agreement is 
found only in the large limit and using a coarse-grained probe. The field theory calculation 
seems to be well-suited for a GGE description. 

Giant gravitons: The half-BPS configurations of A/" = 4 SYM on x R are described by a 
complex matrix QM in a harmonic oscillator potential which can be described by free fermions 
[271 128] in a manner similar to that described in this paper. The gravity duals are the well- 
known LLM geometries [29], which can be described in terms of giant gravitons [30l [311 132] in 
the probe approximation. 

Higher spin duals of 2D CFT's: Recently, [33] has proposed a duality between the large N 
limit of certain 2D CFT's called the minimal models (represented by the coset WZW models 
SU{N)k X SU{N)i/ SU{N)j:_^.i) and a one-parameter higher spin (Vasiliev) theory parametrized 
by A = N/{k + N)fLxed, k, N — )■ oo (see, e.g. the review [31]). The A = limit corresponds 
to 2D relativistic free (complex) fermions and is very similar to the system described in this 
paper. 

D2 branes: The gauge theory system ( |T5l) was considered in [13], where a gravity dual is 
discussed (see Section 5 of that paper) in terms of the near horizon geometry of D2 branes 
compactified on a 'Scherk-Schwarz circle' (i.e. a small spatial circle of size L2 with antiperiodic 
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fermions). The other two directions are periodic, with sizes P,L, which characterize the phase 
diagram of the model is given in Fig. [8] (a). As discussed in [13], the low temperature phase 
P > Per of the model is described by a soliton solution wrapping the temporal circle uniformly, 
while the high temperature ('deconfinement') phase, /3 < Pc, is given by a localized soliton. We 
clearly find evidence of equilibration in the low temperature phase of the field theory in Fig. [9] 
(a). It is also known that the deconfinement phase shows dissipation. It would be interesting 
to understand such equilibration in terms of the solitonic geometries which obviously do not 
have a horizon. 



A/" = 4 super Yang Mills theory: This is clearly one of the most important examples of 
integrable models in AdS/CFT [35]. The giant gravitons described above describe a special 
sector of this model. More generally, operators of A/" = 4 super Yang Mills theory whose 
energies do not scale with A^, can be mapped to states in integrable spin chains which possess 
an infinite number of conserved charges. For energies E ~ 0{N'^) the integrability does not 
persist; thus, a generic perturbation equilibrates, and is described in the gravity dual by black 
hole formation (see, e.g., [M])- The present discussion opens the possibility of equilibration of 
certain perturbations in the integrable subsector and emergence of a GGE; it would clearly be 
of interest to find a gravitational dual of such an ensemble, if any. 

4.3 Black holes and 0{N) entropy 

It is tempting to speculate that the equihbration, discussed above for the matrix QM model, 
and possibly generalizable to other integrable field theories, is described in terms of formation 
of black holes in a dual theory. In many known examples of AdS/CFT, equilibration does have 
a natural interpretation in terms of formation of black holes [^. 

Such a scenario, however, must pass the following stringent tests: 

(i) Black holes are typically not associated with integrable models; while the latter have an 
infinite number of conserved charges, black holes typically have only a few (c/. uniqueness 
theorems). Since a GGE fl37j) is characterized by an infinite number of chemical potentials, any 
black hole describing a GGE must also have an infinite number of parameters. 

(ii) The entropy of the GGE must match the entropy of the known black holes. 

(iii) Besides the static black hole solution, there must be a dynamical collapse solution which 
would capture the equilibration process. The black hole formation time should then match the 
time of relaxation to equilibrium described by the GGE. 

(iv) A non-extremal black hole will Hawking radiate; there must be an interpretation of the 



^''or, in the context of a probe approximation, in terms of the damping of quasinormal fluctuations in a black 
hole background. 
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two-stage process — formation and evaporation — in the field theory description. 

(v) In case we seek a description of equihbration and dissipation in terms of absorption of energy 

by a fixed black hole background in the dual geometry. 

2D: As mentioned above, the c = 1 matrix model ( 12^ . is dual to the 2D bosonic string theory, 
in the linear dilaton background. This theory, as described by its massless fields (the metric, 
dilaton and the massless 'tachyon'), does have the MSWW black hole solution |38l |39] (see 
IIOI HI] for superstring generalizations). Let us discuss the above points to see whether this 
black hole can be relevant for a description of the GGE found in the present paper. 
Point (i): the MSWW black hole solution, is indeed the unique solution of the effective field 
theory of massless modes of 2D string theory, and is specified by only the mass parameter 
However, this solution can be extended to an infinite-parameter generalization in 2D string 
field theory [13]. Each such black hole possesses an infinite number of specified Woo charges 
(corresponding to charges of higher mass non-propagating gauge fields), and could potentially 
correspond to a given GGE with fixed values of the chemical potentials. 

Point (ii): the entropy of the MSWW black hole was found to be (see |12]) of the order of 
^/dstr ~ whereas the entropy of the GGE in Section [2l2l turns out to be of order A^. 
Point (iii): [H] has argued that if we wish to form the MSWW black hole (or its type OB 
counterpart) by collapsing a tachyon shell, physics becomes strongly coupled even before a 
horizon is formed. In the type OB black hole, one can in stead consider collapse of an axion 
shell; in this case, there is no strong coupling problem but the shell is scattered back to infinity 
from near the incipient horizon and a black hole is not formed. 

Point (iv): Reference [H] performs a matrix model computation to rule out any Hawking 
radiation from a possible gravitational collapse of axionic matter. 

Point (v): a fixed black hole background in type OA theory was argued in [ID] to be dual to a 
deformed matrix model [15]. One could in principle consider fermionizing such a model (as in 
( JTOl) ). explore dissipation in this, and attempt to explain this in terms of black hole absorption. 
However, doubts have been raised in [16] about the duality, as the entropy of the matrix model 
is computed there to be 0(1) whereas the entropy of the black hole is 0{N'^). The suggested 
dual matrix model, according to [16], which reproduces an Oi^N"^) entropy, is the one proposed 
in [1?]; however, in this model, the off-diagonal matrix elements play a crucial role, and the 
model is not (at least in any obvious sense) integrable any more. One could also consider the 
proposal in [20] for construction of a black hole background in terms of an asymmetrically filled 
Fermi sea. We leave further studies of these issues for the future. 

^^Thc charged cousins described in [351130] require fluxes which can come from superstring or type theories 
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Higher spin black holes in AdSs: We have already discussed above the duality between 
large 2D CFT's {Wn minimal models) and higher spin (Vasiliev) theories in AdSs. In the 
limit of A = 0, the field theory is that of N number of complex free fermions, and is likely to 
equilibrate as in this paper, with emergence of GGE. As in case of the 2D gravity duals, the 
AdSs also has black holes, the most well-known being the BTZ black hole. Indeed, it has been 
found in |38] that the BTZ black hole has spin three and spin four generalizations. It would 
be very interesting to investigate if this result generalizes to arbitrarily high spin, in which case 
they would be natural candidates for a dual description of GGE. These black holes have an 
entropy of 0{N), matching the central charge c ~ A of the field theory (see [37], e.g.), which 
agree with the 0(A) entropy of the GGE computed in this paper. 

The superstar geometry: In the context of the matrix model-LLM-giant graviton duality 
mentioned earlier, thinned out phase space configurations such as fl72|) in which the fermion 
occupation numbers are fractional have been studied. The resulting LLM geometries are singular 
and are called superstar geometries. These geometries have zero horizon area, but a Bekenstein 
entropy can be assigned to a stretched horizon [15], which turns out to be 0(A) just like the 
standard Boltzman formula for entropy in these systems, and agrees with the GGE estimate of 
entropy. 

D1-D5: The possible relevance of GGE to the two-charge D1-D5 black hole is already men- 
tioned in the previous subsection. The entropy of the two-charge black hole S ~ VQ1Q5 ~ 1/ gstr 
appears to agree with the 0(A) entropy of GGE. 

4.4 Nature of equilibrating observables 

In the discussion of equilibration in an integrable system and GGE, a crucial question (still 
open) is to figure out the nature of observables which indeed equilibrate. E.g. clearly if one 
looks the conserved charges themselves, they do not evolve in time at all. In case of some 
statistical mechanics models, the choice of the variable of interest is dictated by experiment. 
Thus, for both the hardcore bosonic lattice and for the transverse field Ising model, the bosonic 
variables are of more immediate interest. In our matrix model examples, the conserved charges 
can be taken to be the fermion occupation numbers A^; these or 'small operators' of the 
fermion theory, e.g. ([7]), which do not equilibrate. On the other hand, the eigenvalue density, 
as described above, has a simple interpretation in terms of the gravity duals, and it is precisely 
this which shows equilibration. It is an interesting question if equilibration is a generic feature 
of gravity variables in these models. 
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4.5 Topology change and the Gregory-Laflamme phase transition 

The uniform and gapped phases of the double trace model ([3]) correspond to the large L and 
small L sides of the phase boundary AB in Fig. [8] (a). These phases, by definition, refer to the 
uniform and clumped eigenvalue distributions of the spatial Wilson line t/, which arises from the 
gauge of theory of D2 branes wrapped on a small Scherk-Schwarz circle, (see page . What 
do these phases mean in gravity? This question was investigated in detail in [THl EDI ED (see 
[53| [M] for other related references). The Euclidean near-horizon description of such a system 
at large and L (sizes of the temporal and spatial circles, respectively) is that of a D2 soliton. 
If we T-dualize this system along the large spatial circle, we get a Dl-brane soliton which is 
uniformly smeared along a small dual circle L' oc 1/L. This is the description of the phase to 
the right of the line AB in Fig. [8] (a); now as we decrease L, the dual circle L' starts to increase 
and the uniformly smeared Dl soliton starts getting stretched and as L' increases beyond a 
certain critical size, it becomes unstable towards breaking and forming a localized soliton. In 
[55], Gregory and Laflamme predicted such a transition as a thermodynamic phase transition. 
Thus, the dynamical gapless — )■ gapped transition in our model translates to a dynamical GL 
transition from a uniform to a localized soliton. In particular, the observed non-existence of 
such a transition in the matrix model, ultimately due to non-splitting of Fermi liquid droplets, 
implies impossibility of a dynamical GL transition from a 1-brane (uniform soliton) to a 0-brane 
(split solution). 

The original context of the GL transition was in the context of a clumping transition of 
a black string into a black hole. The question about whether this can happen dynamically is 
related to the issue of cosmic censorship and the appearance of a naked singularity, and has 
been discussed extensively in the literature |56l ETJ [Ml EHl EOl Elj. In a supersymmetric variant 
of (fTSll . the phase transition discussed in the previous paragraph would correspond to such a 
transition from a black string to a black hole. Since the technique developed in this paper 
appears to be fairly robust, dynamical transition in this model could share similar features, 
and could shed light on the issue of dynamical appearance of a naked singularity. See also the 
related discussion in page [T2l 
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A Fermionic formulation of MQM 

In this appendix, we will discuss the fermionic formulation of MQM. For a general review, see, 
e.g. [T8l[62]|22 

A.l Single trace model 

Let us first consider the single-trace model ([2]). The singlet sector of this model (see footnote 
[2]) is described by the following A^-fermion action and hamiltonian 

r ^ 

S= dt M--acosft 

^ fc2 

H=J2hi0i,deJ, h = -—dl + V{e), V{e) = acose, h=l/N. (18) 

i=l 

The positions 6i are defined on a circle, and arise from the eigenvalues of U = V diag[e*^']V^~^; 
the fermionic nature arises from the measure (after integrating out the Vi/s) which vanishes for 
coincident eigenvalues. The identification of ^ as 1/A^ comes from the partition function in ([T]). 
A second-quantized version of the above action is written down: 

s= [ dtde i)\e, t) {-ihdt - h{9, de)] ij{e, t), I de ij\e, t)ij{e, t) = n. (19) 



It is useful to describe states of the system in terms of expectation value of the phase space 
(Wigner) distribution 



,p,t) = ldv^^{e + p))^(e-p)e'P^/\ (20) 



We define the position space density as 

N 



i=l 



22 



For specific details of the material presented here, see, e.g. [531 IMl 113 [551 HZl 1551 155] . 
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The normalization condition in implies that 



/ 



d9 p{9,t) = 



dO dp 



■u{9,p,t) = 1. 



(22) 



Large- limit Since h = 1/N, the large limit reduces to 



oo 



Nh= 1. 



(23) 



In this limit, the fermion configurations are represented by a semi classical Fermi liquid with a 
smooth phase space density and position space density; the expectation value u{9,p,t) is either 
or 1 depending on whether the phase space cell (of vanishing area 27rh) is occupied by a 
fermion or not). The regions with u = 1 are called droplets of fermi liquid. 

Phases: Classically, due to the potential V{6) = a cos 6 in (fT8|) the particles tend to clump 
around 6 = tt; quantum mechanically however the states of the system are quantized and Pauli 
exclusion tends to spread the fermions up. The competing tendencies dictate the qualitative 
nature of the ground state: (a) For a > Oc = 7r^/64, the potential is strong and the eigenvalues 
are clumped around 6 = n. Pauli exclusion forces the fermions to occupy levels, but the 
Fermi level is below the top of the potential. This leads to a gap in p{6) around 6 = 0. This 
is called the gapped phase, (b) For a < ac the potential is weak and the eigenvalues are spread 
out (non-uniformly) . The Fermi level is above the top of the potential and there are no gaps in 
p{6). This phase is called a gapless phase. The phase transition at a = is characterized by 
the closing of a gap, when the Fermi level comes up and touches the top of the potential. This 
is a third order phase transition and was found in [2] . 

Double scaling: It is clear from the above that the position of the fermi level p, measured 
from the top of the potential, is 0{a — ac). It turns out that the singularity of the free energy 
is the form Fsing{pN)(x fo{pN)^ + fi + /2(/^^)~^ + .... The double scaling limit is defined by 
A^ — !■ oo, /i — )■ 0, such that pN is held fixed. Our result about the entropy S ~ 0{N) is without 
taking into account such a double scaling. 

The hermitian model: The double scaling limit captures the universal part of the potential 
near the top: 6^0, whose essential features can be equally simply reproduced, by replacing 
cos 6' by (1 — 6'^/2) in ( IT8l) : dropping the 1, the potential becomes V{6) = —a6'^/2, which is 
equivalent to the (singlet sector of the) hermitian matrix model 
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This model is called the c = 1 matrix modelo The double-scaled c = 1 model is mapped 
to two-dimensional bosonic string theory (see the review [18]) by identifying = 1/gstr- A 
similar duality to two-dimensional type theories also exists (see the review |62]). 



Classical droplet motion: Suppose that the Fermi liquid configuration is that of a single 
droplet (e.g. the Fermi sea or a continuous distortion of it), characterized by specifying two 
functions p±{6, t), which demarcate the upper/lower boundaries of the droplet at a specific value 
of 6 at time t, i.e 



(25) 



Using this, we get the following formulae for density, momentum density and energy density: 

p{o,t) = J ^u{e,p,t) = ^ip^{d,t) -p4e,t)) , 

pi9,t)Vi9,t) = J ^pui9,p,t) = ^ipli9,t)-pU0,t))/2, 



H 



d9 dp 



2ii 



h{p,9)u{9,p,t) 



d9 



P-(e,t) 



2ti 



{p'/2 + V{9)), 



d9 



1 -,2 , Tt'^s 



i^pV' + -p' + p{9)V{9) 



(26) 



Here V{9) = a cos 6' for the unitary matrix model ([2]), (fT9|) . whereas V^(^) = —a9^/2 for the her- 
mitian matrix model fl24|) . It turns out that p{9,t),V{9,t) are canonically conjugate variables, 
with Poisson bracket 

{pi9,t),P{9',t)} = de5i9-9'). 

This can be derived from the anticommutation relation of the fermion field, or from the sym- 
plectic Woo structure of the u(6',p)-fields (see [681 169]). The equation of motion of the 'collective 
fields' p, V or p±. follow from the above hamiltonian and Poisson brackets (see comments below 
(pT!) for limitations of the collective field description): 



p{9, t) = de{p{9, t)V{9, t)), V{9, t) = -de 



\V' + \p' + V{9) 



(27) 



These equations of motion (EOM) can alternatively be derived from that of the phase space 
density 



dh 

dt - g^dp + ] u{9,p,t) = 0, 



dp 



(2^ 



^■^Thc periodic boundaries 9 = zLtt need to be supplanted hy 9 = ±9m to make the model well-defined; the 
choice of 9^ is not particularly relevant and does not affect critical/universal quantities. 
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which, viewed from the particle-fixed frame, is simply the particle EOM, and can be regarded 
as the Euler equation or the dissipationless Boltzman equation. This viewpoint provides an 
explicit solution of the equation of motion: 

u{e,p,t) = uoieo{e,p,t),po{e,p,t)), (29) 

where (6*0,^0) is the unique initial phase space point at t = which reaches {6,p) at time t 
through the dynamical evolution ( |T8l) . and uq is the phase space density at t = 0. For the 
hermitian model ( IMll . V = —a6'^/2, and (l29l) becomes 



■>,p,t) = uq ( 6 cosh.(y/at) '^sinh(^/at), —6^/asinh.{^/at) + p cosh.{y/at) ] , (30) 



which can be easily verified to be a solution of f l28|) for h = p'^/2 — a6^/2. In the special case of 
a = (free fermions with zero potential), we get 

uie,p,t) = uo{e-pt,p). (31) 

Note that the collective field equations of motion fl27|) depend on the assumption of a quadratic 
profile ( 125|) . which are clearly violated by 'folds' e.g. in Fig. |6]of Fig. [131 As a result, they are 
not valid for long times (see [681 EH]). In our actual analysis where we are specially interested 
in long time evolutions, we do not use ( 1271) but use ( l28l) or equivalently the particle equations 
of motion in Appendix lEl 

Given the solution f l29l) . and the relation fl2Tl) between p and u, we get the following result 
for p{9,t) at iV cx): 

p{o,t) = J ^uo{eo{e,p,t),po{e,p,t)). (32) 

This can be numerically implemented by sprinkling fermions uniformly in the support of the 
original function uq and let each fermion evolve according to the EOM of ( HM . 

A. 2 Double trace model 

Using the fermionization techniques discussed in the previous subsection, the action (E]) can be 
written as (with h = 1/N, as before) 



S= dt 



2 N 



(33) 

H,nt = j^jdtde de' [i,\e,t)'4j{e,t) cos{e - e')^\e',t)^{e',t)] . (34) 
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The definition of u{6,p,t), p{6,t), V{6,t),p±{6,t) all remain as before. 



H= de 



1 o 



1 



de de'p{e) cos{e- 6') p{e'). (35) 



Because of the limitation of the collective field description in the presence of folds (see comments 
below fl3T]) ) for computation of the time evolution we use the equation of motion of the phase 
space density or equivalently the particle equations of motion as in Appendix [El 

B Integrability and GGE 

The single-trace model (1181) . (fT9|) is clearly integrable, since it consists of free particles. 
Classically one can consider the energy ej of each particle, ^ = 1, A^, to be an 'action' variable. 
Quantum mechanically, each linearly independent A^-fermion state can be specified by saying 
which single-particle energy levels, e^, are occupied. The quantities 

N 

1=1 

are all constant in time, and commute with each other. For finite A^, only A^ of them, Ip 
{p = 1, A^) is enough to determine the energy levels (up to some discrete ambiguities which 
we shall ignore), hence In+i, In+2, ■■■ are all dependent on the first A^ Jp's. In the N ^ oo limit, 
of course, the charges ( 136|) are all independent. 

We can define the Generalized Gibbs Ensemble (GGE) in terms of the A^ conserved charges 
as follows: 

^•gge = -77^ exp [- Ppip] , Zqge = Tr exp [- V f3plp] . (37) 

^GGE 

P P 

For more details, see [1]. 

In the present discussion, it is convenient to introduce the fermion occupation numbers Nm, 
m = 1, oo, where Nm = 1 if the single particle energy level is occupied, and = if it is 
not (note that our fermions are spinless, and that the energy levels are non-degenerate). These 
can be defined in terms of the creation/annnihilation operators through the equations 

iPie) = ^c„0™(e), h{e,de)<Pm{G) = emM^), (38) 

m 

Nm = clCm- (39) 

In terms of the A^'s, Ip = Ylim^m^m- For A^ = oo, in stead of ( 136|) . one can use A'm's as the 
(independent) conserved charges, which leads to the following definition of the GGE: 

^•gge = exp[- ^ /i^A"™], ^gge = Tr exp[- ^ p-mNm]- (40) 



GGE 
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With this definition, we get, by exphcit evaluation (see [70]). 

InZ = Vln(l + 6-'^-), (iV„)GGE = -^InZ = — i— , (41) 

m ^ 

which gives a remarkably simple Fermi-Dirac distribution which involves only one chemical 
potential. The chemical potentials are now trivially found: 

f^m = In (— i l) . (42) 

The entropy of the GGE is given by 

oo 

'S'gGE = -Tr^GGE In ^iGGE = -J]] [(^m)GGEln(iVm)GGE + (1 " (^''m)GGE) ln(l - (iV;m)GGE)] • 

m=l 

(43) 

The GGE hypothesis is that, if we start from an appropriate excited state |\l/(0)), there exists 
a certain class of observables O such that, as t — > oo 

(vl/(t)|0|vl/(t))-^Tr(f?GGEO), (44) 

where the GGE is defined by the charges 

(iV™)GGE = (*(0)|iV^|*(0)) = {^{t)\Nm\^{t)), (45) 

or equivalently by the chemical potentials fim, determined through (H2|) . 

Note that if the initial state |^I^(0)) is an eigenstate of the number operators Nm, then by 
(H5l) . (Nm) = or 1. The GGE constructed from such an initial state has vanishing entropy 
'S'gge = 0, according to ( l43l) . 

B.l Verification of the GGE hypothesis for p(9) 
We wish to verify (jH]) where the observable is the density variable 

o{t) = p{e,t) = 5^4(t)c.(t)C(e)0„(0), 

where Cm(t) = exp[iem,V^]; = c]„exp[— iem^/^]- The LHS of (l44l) becomes, in the 

Heisenberg picture, 

{p{e,t)) = ^{cl,Cn)(j)l,{x)(j)n{x) exp[-i{em - en)t/h], (46) 
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where (...) = (\E'(0)|...|\l/(0)) represents expectation value in the initial state |\E'(0)). We will 
evaluate this expression in (156|) where the initial state is prepared according to QQD. The large 
N limit of f H6|) is given by the semi classical expression fl32|) . 
To compute the RHS of ( 144|) for O = p{0), first note that 

GGE ^mnij^rn) GGE, (47) 

where {...)gge = Tr(^GGE[---])- Using this and fl38|) . we get 

(p(^))gGE = J2^Nr,,)GGE\4>.^m'■ (48) 
m 

Note that this is independent of time, as expected. Here {Nm)GGE are to be determined from 
their values in the initial state |^E'(0)) 

(iV„^)GGE = (^(0)|iV™|*(0)). (49) 

Thus, 

(p(^))gge = J](v^(O)|iV^|v^(O))|0„,(^)|^ (50) 

m 

Since these uniquely determine the chemical potentials through ( H2|) . the GGE is specified 
completely. 

Verifying the GGE hypothesis amounts to showing that (l46ll tends to ( l50l) as t — t- oo. In 
Section [C] we will compute both these expressions, and in Section 12.11 present the result in 
support of the hypothesis. 

B.2 Phase space density in GGE 

In this section, we will compute {u{0,p))gge, where u{0,p) is the Wigner phase space density 
(|20|) . By using (!38|) and ( 147|) . it is easy to show that 

{u{9,P))gGE = '^{Nm)GGEUmid,p), 
m 

uUO, P) = Jdrj C {O + I) 0^ {9 - I) e^^^/^ (51) 

where (A^m)GGE is given by f HOj) . In Section [DTTI we evaluate this expression in the semiclassical 
limit. 
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C Quantum quench dynamics (QQD) 



Suppose we consider a sudden change of the parameter a in the single-particle hamiltonian f[T^ . 
from at t < to a/- at t > 0. Time evolution works as follows. We consider the ground state 
\F') of the system (Fermi sea in this case) appropriate to the parameter for t < 0. Quantum 
quench dynamics refers to the sudden approximation in which the system evolves for t > 
according to the new parameter aj, from an initial state equal to the old Fermi sea 

1^(0)) = \F'). (52) 

Clearly the system is integrable after t = and our discussion in the previous section applies 
(we need to identify the a parameter with a/). 
We wish to compute the time-dependence of 

{p{e,t)) ^ {F'\p{e,t)\F') = {mmim), m 

for t > 0. Here the second expression is in the Heisenberg picture (in which the state remains 
the initial state (152|) ). while the third expression is in the Schrodinger picture in which |^I/(t)) 
is the time evolution of (132]) according to the new Hamiltonian. 

To proceed, let us express the old Fermi sea \F') in terms of eigenstates in terms of the new 
number operators Nm- Suppose we denote the eigensystem of the old hamiltonian (with a^) as 
e^. The functions 0^ can clearly be written in the new eiganbasis thus 



m m 

Just like (l38l) . the Fermi field can be equally well expanded in terms of 0^; thus 

n m 

we get 

— B^^Cn, — -BmnCli, Cn = C^Bmm = C^B^^. (55) 



Using the expression (H6l) we get 

oo N 

{p{e,t))= J2 J2^rmBrn<l>*Jx)U^)exp[-t{e^-e^)t/h], (56) 

m,n=l r=l 

where we have used ( 155|) and the defining property of the old Fermi-sea: 

(F'l^VjF') = iV;„ K, = fJnk, /„ = l(n = l,...,iV), /„ = 0(n>iV). (57) 



We explore ( 1561) numerically in the text for various values of A^. The — ?■ oo limit can be 
computed in the semiclassically in the droplet picture as in fl32l) . 
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C.l GGE for QQD 

We wish to compute (p(^))gge in the context of QQD, using the general expression ( 1501) . for 
the initial state fl52|) . 

The expression (149|) becomes 

N 

{NJgge = {F'\clcm\F') = {F'\c^B:y,Bkm\F') = (58) 

n=l 

where we have used (157|) and (!55|) . Using this expression in (148|) . we now get 

AT 

5^5^|5„„,|V™Wr (59) 

m n=l 

In Section 1271] we verify the hypothesis (144|) for the fermion density p{0,t) (more precisely, 
its various moments (jl])), by comparing the large time behaviour of fl56|) (or f l32|) ) and fl59|) . 

C.2 Adiabatic time evolution 

Suppose we have a time-dependent Hamiltonian H(t), e.g. ( fT9l) with a = a(t)-parameter, or 
(IMl) with ^ = ,^(t)-parameter. Suppose the time-dependence starts after some time to; we wish 
to study the time-evolution of a wavefunction |^I^(t)), which at is an eigenstate: iJ(to)|^(to)) 
= i?o|^(to))- QQD is an example in which the time- variation is in the form a step function at 
t = 0. An adiabatic time evolution, on the other hand, refers to a sufficiently slow time- variation 
such that one can approximate {"^(t)) as the instantaneous eigenstate of H{t) with eigenvalue 
E{t); this notion is well-defined in the absence of level crossing, which can be ensured if the 
spectrum is discrete which sets a finite time scale. In the examples considered in this paper, we 
are often interested in the time-evolution of the ground state, and by an adiabatic expectation 
value Vadiabit) we will always mean expectation value in the instantaneous ground state finst(^) 

ground,inst gTound,instify ) \. I 

In the context of ( IT9|) . the ground state of a specific value a{ti) = ai is clearly given by, 
in an obvious notation where we denote all instantaneous quantities by an additional subscript 
(•••)i 

N 

(-^l|Cl,mCl,n|-^l) = Smnfl,m, fl,m = ^nm- (60) 



n=l 



By expanding ip{6) in this new basis, we get 



N 



{P{e))ad.a, ^ {Fl,\pmK) = E \^'mm'- (61) 



m — 1 



^'*At a critical point, the spectrum is gapless; there we will simply define Vadiab{t) = Winst(0- 
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Figure 13: Equilibration of the phase space density in the V{U) = case. The perturbation (j62|) with 71 = 1 
at t = (the left figure) evolves to the right configuration. 



The time dependence of the adiabatic quantity is reflected in the appearance of the instantaneous 
eigenf unctions. 

D An analytic calculation of relaxation 

In this section, we will consider time evolution in the single-trace model ([2]) with a = 0. (This 
model is also equivalent to = in (I3l).) In this case, the ground state is given by p±{0) = ±1/2 
(see (l26l) for notation) @ Suppose we perturb the ground state configuration at t = by a 
sinusoidal deformation 

p+{e,t = 0) = ^ + bcosinO), {\b\ < 1/2), p-{0,t = 0) = -p+{e,t = 0), (62) 

where b is the amplitude of the perturbation (see Fig. [T3|) . We wish to evaluate the time 
evolution of this perturbation using the action ( IT9|) with a = 0. 

We will evolve the configuration as in ( l3T]) . which essentially says that ^4(6*0 + Pot,po) = 
Uo{9q,pq)] in other words, to get any point of the evolved Fermi surface at time t, we should 
simply evolve points on the initial Fermi surface according to the equation of motion of ( !T8|) 
with a = 0, which is given by 

p = 0, 9 = p. (63) 

In order to solve this equation for points on the Fermi surface, it is convenient to employ a 
parameter s E [—n,n] and rewrite the initial configuration (l62l) as0 

p+(s,t = 0) = ^ + 6cos(ns), e{s,t = 0) = s. (64) 



^^Note that this corresponds to p{6) — ^,7^(0) = 0, which obviously satisfies (P7)) for V = 0. 
^^From now, we omit since we can derive it trivially from p+. 
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Then from fl63p we obtain the Fermi surface at time t as 

t) = - + 6cos(ns), 9{s,t) = s + [ - + bcos{ns) \ t. (65) 



This solution with n = 1 is plotted in the middle of Fig. [131 the 'folds' appear because of the 
stretching caused by different speeds of fermions at different heights from the ^-axis. As time 
progresses, more or more 'folds' keep appearing. 

In the perturbation (|62l) . only p„ is excited and other pm are zero. We evaluate how these 
modes behave as t increases. From the solution fl65|) . we calculate the m-th moment as 



r de , , r dsde(s,t) ,,,,,, 

Pm{t)=J — cos(m0) (p+ = y — cos(mg(g,t)) (p+ 

= J ^ (1 — nfttsin (ras)) cos ^ms + mt ^- + &cos(ns)^^ ^- + bcos{ns] 

= — J (iscos(ms) cos ^mt ^- + 6cos(?T,s)^^ ^- + 6cos(?7,s)^ 

H — J ds nbtsm{ns) sin(ms) sin ^mt ^- + 6cos(?7,s)^^ ^- + 6cos(ns)^ . (66) 

These integrals will be given by the Bessel function J^. For example, if n = 1, the integral ( 166|) 
with m = 1 and 2 are evaluated as 

pi(t) = jMbt) cos j , p2{t) = -jJ2{2bt) sin (t) . (67) 



For large t, they behave as 

■ ^3/2 Y ^5^^^ \^"'' 4 y -^'^ \^4^ 

and 



Pi(^) ~^ "T^os ( ^'^ ~ ~r ) ( 1~ ) ' (^ oo), (68) 



/^2W ^ ^^cos(26t-^)sin(t), (t^oo). (69) 

Other Pm{t) also exhibits a similar decay. In Section [273] we will compare the time-evolution of 
various modes and find a picture of an energy cascade in which lower frequency modes die out 
before the higher frequency modes, causing a transfer of energy from higher to lower frequencies, 
in a manner similar to energy cascades in turbulence. 



^''The expression in the first line in this equation is crude, since the droplet has many folds at a late time as 
shown in the center of Fig. [iSjand we need to subtract the white region below p+. However after the change of 
the variable from 9 to s, this subtraction is automatically involved. 
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D.l Asymptotic and GGE form of the Fermi surface 



As we mentioned below (l65l) . the asymptotic form of the droplet consists of an infinite number 
of 'folds' which eventually densely fill an entire band. When suitably coarse-grained (the 'grain' 
size can be finer and finer as time progresses), this describes a thinned out Fermi liquid, as 
schematically plotted on the right of Fig. [T31 

Since we have already verified the GGE hypothesis of our system for the density variable, 
let us describe the above 'thinned out' Fermi liquid in terms of the GGE. We begin with the 
GGE value (|5T1) of the phase space density. The single-particle Hilbert space can be labelled by 
the momentum wavefunctions 

which trivially gives Um{0,p) = S{m — p/h), and hence by ( l5Ti) (keeping the notation (...)gge 
implicit) 

ui9,p) = y^N^5{m-p/h) = N,,n = (70) 
^ — ' Zn 

m 

We can identify p{p) as the momentum-space density 

p{p) = j d9u{9,p). (71) 

We can compute AT^'s from a given initial quantum state by using (149|) . However, the calculation 
is much better done semiclassically, when the initial state is given in terms of a classical droplet 
configuration ( l62l) . Note, first of all, that the momentum-space density ( 17T|) in the present case 
is a constant of motion. This is easy to see since h = pV2 and all fermions have a constant 
momentum so the momentum density remains constanto We can, therefore, compute the 
momentum density from the initial configuration ( 162|) . Note that from ( l62l) and ( 125|) we can 
find out the initial value of u{9,p) and hence p{p) by using (ITTil . The result of this calculation 
is 
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1 



u{9,p) = —p{p)= ^1 bl<i-&, (72) 

^ cos""*^ (p — I)) otherwise, 

which gives us a quantitative representation of a 'thinned out' Fermi liquid on the right panel 
of Fig. 



28 It can also be seen from (gS]) and ((7T|) . 

2^ Geometrically the RHS of ([7T|) is given by d+{p) — 0~{p) where 0±{p) are the counterparts of p±{0) in (pS)) . 
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E Details of our numerical analysis 



We explain the details of our computations in Section |2] and [3l 

Exact computation We use an exact computation for the single trace matrix model ([2]) at 
finite A^. In this case, the eigen function (pmid) is given by the Mathieu function. We numerically 
evaluate this function by using GCC with the GNU Scientific Library. 

Semiclassical computation For the single trace matrix model (j2]) at large- and the dou- 
ble trace matrix model (jS]), we use a semiclassical approximation. At the strict large- A^, the 
dynamics of the matrix models are described by the droplet u{6,p,t) in the phase space. We 
discretize the support of the initial droplet by N uniformly distributed points (we take N to be 
0(10'^).) Each of these points is then evolved according to the classical equation of motion 

Oi = Pi, Pi = asm6i, (73) 

for the single trace model ([2]), and according to the classical equation of motion 

0^ = Pi, P^ = -^Y1 ~ ^j)' ('^^) 

j 

for the double trace model ([3]). Here 9i and pi {i = 1, ■ ■ ■ , N) are the points of the discretized 
droplet. Evolving these points in this fashion is equivalent to a discretization of the Euler 
equation ([281) . By using the Runge-Kutta method, we solve these equations with a given initial 
conditionlfj. Note that the above use of the classical equations means that we are in the strict 
large- A^ limit. As a result, the number N, which is simply a measure of discretization of the 
classical phase space, is not related to A^ in the original matrix model (while 1/A^ measures 
quantum effects, 1/A^ measures some discretization error of classical equations). 
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